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Superconducting, flux-based qubits are promising candidates for the construction of a large scale quantum 
computer. We present an explicit quantum mechanical calculation of the coherent behavior of a flux based 
quantum bit in a noisy experimental environment such as an integrated circuit containing bias and control elec- 
tronics. We show that non-thermal noise sources, such as bias current fluctuations and magnetic coupling to 
nearby active control circuits, will cause decoherence of a flux-based qubit on a timescale comparable to recent 
experimental coherence time measurements. 
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L INTRODUCTION 

The possibility of observing quantum coherent behavior in 
macroscopic devices, such as superconducting circuits, was 
first suggested by A. J. Leggett in the mid 1980's [1]. Since 
then, experimentalists have attempted to observe such behav- 
ior Q|3|| and, more recently, to extend the idea to a full scale 
quantum computer capable of executing quantum algorithms 
such as Shor's algorithm [4]. Superconducting qubit types in- 
clude charge qubits: created by a superposition of the number 
of electrons on a superconducting island |5]; phase qubits: the 
superposition of energy states in a single Josephson junction 
@; and flux qubits: the superposition of the quantity of flux 
threading a superconducting ring interrupted by a Josephson 
junction Q. Quantum mechanical behavior has been veri- 
fied in all three of these qubit systems, and experiments have 
shown that each system can be prepared in quantum super- 
positions. The coherence times which result from these en- 
vironmental interactions, however, are far shorter than those 
predicted by current theories of open quantum systems. 

Non-thermal sources of noise, such as fluctuations from 
room temperature laboratory control and measurement equip- 
ment and active circuit elements operating in the vicinity of 
the quantum device will contribute to the decoherence of the 
qubit in ways that are unique to different experimental con- 
figurations and methods. In this work we calculate the co- 
herence time of a flux-based quantum bit exposed to classi- 
cal noise sources consistent with those that would be found 
in an integrated circuit environment. We explicitly calculate 
the evolution of the qubit wavefunction under the influence of 
a Hamiltonian including environmental noise and infer deco- 
herence times via an ensemble average of such calculations. 

The future of solid state quantum computation will likely 
involve the integration of quantum bits onto a monolithic cir- 
cuit also containing classical electronics used for quantum 
state preparation, manipulation and readout |4, 8]. For a large- 
scale quantum computer the classical control electronics will 
most likely be digital and be implemented in a well under- 
stood, robust technology such as superconducting rapid sin- 
gle flux quantum (RSFQ) technology. RSFQ logic is an in- 
tegrated circuit family that uses single flux quantum (SFQ) 
magnetic pulses as data bits. A diagram of a likely physical 
layout of an integrated circuit is illustrated in Fig. [2 where 
RSFQ digital circuitry is placed near by and inductively cou- 
pled to a rf-SQUID qubit. Unlike charge and phase based 
Josephson qubits, the flux based rf-SQUID qubit is induc- 
tively coupled to the environment making it susceptible to the 




FIG. 1: A RSFQ circuit (left) physically near and inductively cou- 
pled to a rf-SQUID qubit (right). The RSFQ circuit passes digital 
bits as quantized units of magnetic flux. The physical proximity of 
the active circuit will destroy the coherent behavior of the rf-SQUID 
qubit. The elliptical lines surrounding the RSFQ inductor are the 
field lines created by a passing data pulse. 



effects of stray magnetic fields. The quantum dynamics of 
the rf-SQUID is exceedingly sensitive to the applied exter- 
nal magnetic field. The mutual inductance of the rf-SQUID 
qubit facilitates simple coupling procedures between qubits 
and classical circuitry, but also leaves the qubit vulnerable to 
unwanted coupling to active circuit elements integrated on the 
same chip. Although one clearly will attempt to shield the 
qubit from such stray fields these measures cannot be perfect. 
Furthermore, in order to be of any practical use, the flux qubits 
can not be completely isolated. Studies have been performed 
evaluating the effect of mutual inductive coupling in standard 
CMOS circuit technology |9]. Due to the robustness of digital 
logic, classical circuitry has a high level of immunity to these 
types of effects, contrary to quantum coherent technologies, 
which will be extremely sensitive to noise coupled through a 
mutual inductance. 

Though quantum dynamics of the rf-SQUID have been ver- 
ified, the macroscopic quantum coherent (MQC) oscillations 
of flux, predicted by A.J. Leggett 1 1] have not yet been di- 
rectly observed due to the lack of sufficiently sensitive and fast 
means of measurement. By integrating the entire MQC exper- 
iment onto a monolithic circuit we believe it will be possible 
to create the necessary electronics to conduct the experiment. 
Furthermore, such chips used to perform the necessary exper- 
imental procedures can be used to explore decoherence mech- 
anisms arising from qubit coupling to classical noise sources. 

The calculation performed in the following sections models 
the decoherence of a rf-SQUID qubit integrated with active, 
classical circuitry. The effect of the classical circuits on the 
qubit is modeled by the inclusion of fluctuating terms in the 
Hamiltonian of the rf-SQUID. The classical sources of noise 
are characterized by Gaussian random variables with a finite 
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bandwidth and predetermined power spectral density (PSD). 
The predicted qubit coherence times depend upon character- 
istics of the classical noise that may be measured experimen- 
tally. 



II. THE DECOHERENCE CALCULATION 

Decoherence can most easily be defined as the deviation 
of the behavior of a quantum mechanical system from that 
predicted by the Schrdinger equation for the closed quantum 
system. Traditionally, decoherence is defined in the context of 
the density matrix (p) representation where loss of coherence 
in a quantum system is indicated by the suppression of the 
off-diagonal elements of the density matrix. In a given basis, 

1 N 

t=l 

where N is the number of individual systems composing the 
ensemble about which quantum statistics are being described 
and \ipi(t)) is the state of the i th system. When the off- 
diagonal elements of p are completely suppressed evaluation 
of the Schrdinger equation for the closed system ceases to 
be an accurate predition of the evolution of the state \ip) of 
the open system. Decoherence models have been created to 
predict the time scales on which non-deterministic degrees of 
freedom of the environment become entangled with the qubit 
thereby destroying the coherence of the quantum system I Kill . 
There is a large body of work on the decoherence of a sgin-1/2 
system coupled to a reservoir of harmonic oscillators | lll llal . 
To our knowledge, however, an explicit, time-dependent cal- 
culation of the evolution of a flux-based qubit including the 
effects of a noisy environment has not been completed. 

Two distinguishable decoherence mechanisms contribute to 
the suppression of the off-diagonal elements of p. Relax- 
ation is associated with an increase or reduction of the ex- 
pectation value of the energy. Dephasing is an adiabatic pro- 
cess whereby the phase of the system wavefunction becomes 
randomized. Dephasing typically occurs on a much shorter 
timescale than relaxation fl3ll . 

The Hamiltonian for the rf-SQUID, derived in 1 1 ], is given 

by, 
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In this equation C is the capacitance of the Josephson junc- 
tion, $3 is the magnetic flux applied to the device externally 
and L is the loop inductance. The independent variable, $ 
is the flux threading the superconducting loop. The first term 
represents the kinetic energy of the SQUID, while the second 
and third terms constitute the potential energy of the super- 
conducting inductor and Josephson junction, respectively. 

The qubit potential has two independent degrees of free- 
dom, the height of the barrier separating the two minima, and 
the relative depth of the two minima. The Hamiltonian can be 
rewritten in a form similar to that of a two-state system as, 
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where A is the tunnelling matrix element, and e is the differ- 
ence in energy between the ground states of the two wells. The 
operators a x and a z are Pauli matrices. In principle these two 
degrees of freedom are separable, and can be considered inde- 
pendently. In a laboratory setting, however, they are subject to 
the same environmental influences and should be considered 
simultaneously when examining decoherence. The flux states 
of the system correspond to the eigenstates of a z . In this basis 
fluctuations in the barrier height are referred to as a x fluctua- 
tions because they modulate the energy level spacing and fluc- 
tuations in the relative depths of the wells are called a z fluc- 
tuations. When the flux bias deviates from $0/2 by greater 
than 10 -4 <l>o, the energy and flux basis states are nearly the 
same. However, when the system is flux biased at exactly 
$o/2 the energy bases are non-local in flux; there is finite 
probability of finding the flux in either well. Traditionally, a z 
fluctuations are considered to be the most destructive to the 
coherence of a system. In a typical experiment[2], the sin- 
gle rf-SQUID junction is replaced by a double junction loop 
of small self-inductance. This added inductance will promote 
coupling between the environmental fluctuations and a x de- 
gree of freedom in the system. The contribution to dephasing 
from a x and a z degrees of freedom will depend on the ratio of 
the self-inductances of the junction loop and SQUID loop, re- 
spectively. Reducing the junction loop inductance will allow 
a x dephasing to be ignored. 

As stated earlier, the calculation performed in this paper is 
an explicit solution of the time-dependant Schrdinger equa- 
tion. Previous calculations of the decoherence of a rf-SQUID 
have focused on modeling the interaction of the environment 
as a continuous weak measurement I14ll or by reducing the 
problem to a two-state system and using the spin-boson for- 
malism 1U2I. IT5ll . The shortcomings of these models are that 
they approximate the effect of a general environment on the 
system, but do not accurately reflect the dominant sources of 
noise in a circuit environment. For example, the spin-boson 
formalism reduces the qubit to a system with only two energy 
levels and models the interaction between the qubit and envi- 
ronment as bilateral, in that the state of the qubit affects the 
state of the environment. In an integrated circuit environment 
the back action of the qubit on the sources of the fluctuating 
fields is insignificant. 

The illustration in Fig. ^depicts an RSFQ integrated cir- 
cuit, consisting of resistors and superconducting inductors and 
Josephson junctions near a rf-SQUID qubit. The magnetic 
data pulses that propagate through the RSFQ circuit will in- 
ductively couple flux to the qubit perturbing the otherwise 
constant external flux bias; it is necessary to determine the 
magnitude of these decoherence causing fluctuations. Stray 
magnetic flux coupled to the qubit from the RSFQ circuitry is 
determined by the magnitude of the current pulses in nearby 
RSFQ circuits and the mutual inductance (M) between the cir- 
cuits and qubit. A typical estimate for M when the circuit and 
qubit are 20/im apart is on the order of IQfH. A SFQ data 
pulse passing through a superconducting inductor of typical 
size in RSFQ circuits has a peak amplitude of approximately 
150/i/l and a duration of 2 — Aps. With these parameters the 
excess flux coupled to the qubit is 720/i<I > o- An analysis of the 
decoherence effects in such a situation follows. 

The effect of a large number of nearby RSFQ circuits, each 
passing SFQ pulses at seemingly random times from the point 
of view of the qubit, may be represented by introducing a ran- 
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dom component to the flux in the vicinity of the qubit. Thus, 
our calculation was performed by introducing noise terms di- 
rectly into the Hamiltonian, Eq. |2 and then solving the time- 
dependant Schrdinger equation with the random Hamiltonian. 
The rf-SQUID can be perturbed by the flux environment in 
two ways: the external flux bias applied to the qubit can de- 
viate from $o/2, or magnetic flux can couple to the Joseph- 
son junction thereby changing the effective I c . As shown ear- 
lier the influence of the fluctuations of the applied flux bias, 
a z noise, dominates the effect of qubit junction critical cur- 
rent fluctuations, the a x noise. Therefore a random term was 
added to in the Hamiltonian and the evolution of the qubit was 
computed. The initial wavefunction, used in the calcu- 
lation was determined by the ground state of the system with 
slightly less than $o/2 applied flux. At milliKelvin temper- 
atures the system will completely relax to this state. |\&o) is 
then projected onto the basis functions of the symmetric dou- 
ble well potential. The state of the system is nearly a pure 
superposition of the ground and first excited state of the sym- 
metric potential. There are, however, excitations of higher 
lying states that are small but nonetheless included in the cal- 
culation. 

The time evolution of the system is determined using a time 
varying propagator operating on the wavefunction, 

\ip m (t + At)) = exp(iH m (t + At)At)|Vm(*)>, (4) 

where H m (t + At) is the Hamiltonian of the system with 
$ x (t + At) and I c (t + At) ; we also chose At to be lOOps, 
much shorter than the period of the MQC oscillations. At each 
time step the wave function \ip m (t)) is projected onto the basis 
states of the new Hamiltonian H m (t + At), which are deter- 
mined before the calculation and stored in a look up table, 
and then evolved through time At to produce the new wave- 
function, \ip m {t + At)). The result for short time-scales is a 
randomization of the phase of the coherent oscillations of the 
flux. On a longer time-scales the populations of the energy 
levels will deviate from those of l^o)- Since the loss of phase 
coherence occurs on a shorter time-scale that is the dominant 
form of decoherence and thus will be the subject of this work. 

In order to calculate the dephasing time of the MQC os- 
cillations in the fluctuating potential the simulation described 
above is repeated many times, and the separate results are av- 
eraged as shown in Eq. ^ An experimental technique simi- 
lar to the one described here was used in 1 16] to measure the 
phase stability of classical oscillators. For our calculation we 
used N = 50. For each trial the phase of the oscillations 
becomes unpredictable after some period of time. Thus av- 
eraging over a large number of random phases leads to the 
vanishing of the off-diagonal elements of the density matrix 
inEq.G] 

A set of N evolutions was performed and summed for each 
value of flux noise variance, and the results are plotted in the 
surface plot of Fig. |2] The bandwidth of the noise in Fig. |2] 
is AGHz. To quantify the dephasing time resulting from each 
value of flux noise variance, the damped sinusoid resulting 
from the set of N evolutions was fit to a function of the form, 

P(t) = exp^D^cos^ot + 7r), (5) 

where ojq is the resonant frequency of the tunnelling flux, 
which is also equal to the energy splitting between the ground 
and first excited states. The value of the dephasing time 




FIG. 2: A surface plot of the data obtained from the calculation of 
decoherence of the rf-SQUID qubit. As the variance of the flux noise 
influencing the qubit is increased, the coherent oscillations of flux 
last for shorter amounts of time. 



constant, D^, for each value of flux noise variance was 
determined using a least squares fit to the data obtained in 
Fig. |2 Fig. |3 shows the plot of as a function of 
the flux noise variance, a^, applied to the qubit. From the 
graph in Fig. [3] the relation between and the variance is 
Z?0 = 3.63 • 1CP 4 (ct^) 2 (1/ns), where er^ is in units of /i$o- 
Since <r? = S^>(u>)duj, D,p is linearly proportional to the 
amplitude of the power spectral density of the environmen- 
tal flux noise at a given bandwidth. Fig. |2 shows that if the 
flux noise variance is increased above approximately 10/i$o 
, then the coherence of the qubit lasts less than 20ns. From 
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FIG. 3: A plot of the dephasing rate as a function of flux noise vari- 
ance. The bandwidth of the noise signal used for this data is AGHz. 
Error bars are provided for selected data points indicating the confi- 
dence of the fit. 
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the physical scenario described above, it is evident that M 
must be reduced by at least two orders of magnitude, to below 
■ lfH, in order for the excess flux to be reduced to these levels 
and for acceptable coherence times to be achieved. 

An analysis was also performed to measure the dependence 
of Z?0 on the bandwidth of the noise signal. It was seen that if 
the noise is bandlimited to a cutoff frequency, uj c , below uio, 
then varies approximately linearly with the noise band- 
width. If uj c > uj , is nearly independent of the bandwidth 
of the noise signal. When the noise is bandlimited so that 
uj c < uiq the frequency of the oscillations is no longer ujo due 
to the fact that on the timescales in question noise bandlimited 
below ujo can no longer be considered white, and Q x 
This yields a noticeable increase in the coherent tunnelling 
frequency of the system. At noise frequencies greater than ui n 
the noise can be considered white and <& x = ^ over the 35 
ns time scale of the simulation. At large noise frequencies the 
dephasing rate, D^, is dependent primarily on the variance 
of the flux noise and nearly independent of the bandwidth of 
the noise. In an integrated circuit environment, especially one 
operating at GHz speeds, the bandwidth of the noise created 
by the circuitry will certainly exceed that of uj , which in this 
case is approximately 280 MHz. 



III. CONCLUSIONS 

The ultimate goal of integrating classical electronics with 
coherent quantum circuitry is to efficiently prepare and ma- 
nipulate the coherent states of the qubits. This will require that 
classical circuitry, such as RSFQ, be placed in close proximity 



to the quantum coherent bits in order to provide the necessary 
interaction, which in this case is inductive. Unwanted interac- 
tions between unintentionally coupled classical circuitry and 
qubits will no doubt exist and it must be determined what ef- 
fect that will have on the coherence time of the qubits. We 
have calculated that for current RSFQ technology there is a 
threshold value for the mutual inductance between RSFQ cir- 
cuitry and a qubit of approximately 0. lfH below which qubit 
coherence times may persist for longer than 20ns. 

We have used an explicit formulation of the Hamiltonian 
of the rf-SQUID in a noisy environment in order to determine 
a realistic value for the dephasing time in an integrated cir- 
cuit environment. Moreover, the parameter that accounts for 
the loss of coherence in this model, the flux noise variance, 
is easily measured by incorporating a dc-SQUID magnetome- 
ter in the vicinity of the qubit. The approach presented here 
for estimating the coherence time of a rf-SQUID qubit is not 
intended as a general model of decoherence for two-state sys- 
tems. Rather it provides an explicit model from which co- 
herence times of rf-SQUID qubits may be estimated based 
upon experimentally measurable quantities. Models such as 
this one coupled with measurements of the flux noise at the 
location of the qubit will help quantum computer architects 
design large scale computers capable of executing the algo- 
rithms for which they were originally intended. 
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